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Our work has argued for a particular scaling form governing the distribution M(£^,t) of magneti- 
sation over bias ^, for a system of dipolar-interacting molecular spins. This form, which was found 
in Monte Carlo (MC) simulations, leads inevitably to a short-time form ~ t^^'^ for the magnetisation 
relaxation in the system. The authors of the Comment argue that the magnetisation should decay 
rather as ~ t"^ , with the exponent p depending on the lattice type- and they argue this form is 
valid up to infinite times. They also claim that our conclusion is based on an assumed exponential 
dependence of the function M{£^,t) on Tde{Oy the effective molecular relaxation time. In fact our 
results do not depend on any such dependence, which was used merely for illustrative purposes, 
but only on the scaling form we found. Repeating our MC simulations for different lattice types 
and different parameters, we always find a square root relaxation for short times. We find that the 
results of the comment are flawed because they try to fit their results over far too large a range of 
times (including the infinite time limit, where no simple theory applies). 
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The comment— misrepresents both our work and the 
physics of the problem. Our results are essentially that 
for both strongly polarised^ and depolarised^ dipolar sys- 
tems, the short-time relaxation (after an initial tran- 
sient) has the form SAI ^ ^/t, regardless of lattice type- 
with concomitant results for the hole width and the scal- 
ing function Af(^,t) introduced in^. It is argued in the 
comment^ that (i) our results are based on the assump- 
tion of an exponential form for the function M(^, t); and 
(ii) that the magnetisation relaxation in lattices of vari- 
ous symmetry may be described over huge time ranges, 
encompassing up to 3 orders of magnitude in magnetisa- 
tion (from Af = 1 to Af = 10""^) by a power law form 
dM{t) ^ tP, with p dependent on lattice symmetry. We 
respond to these in turn: 

(i) The effective relaxation rate tJ^^(C) was introduced 
in^, and extracted and studied using direct i-dependent 
MC simulations. It describes molecules in large bias 
fields ^, and, we claimed, has a Lorentzian-tail depen- 
dence tJ^^(^) ^ As before we write the total nor- 
malised magnetisation of the system M{t) — J d^AI{£^, t), 
where M{^,t) is the normalised distribution of magneti- 
sation over the local bias fields ^ acting on each molecu- 
lar spin; in terms of the probability Pa-{£,,t) for a spin 
to be in a bias field ^ at time t, one has M{^, t) = 
[P|(^, t) — Px(^, t)]. In commenting on our paper Alonso 
and Fernandez have introduced a function f{^,t) which 
is nothing but f{£,,t) — —M{^,t). They argue that 
our conclusions are based on the assumption that our 
M{^,t) depends on Tde{£.) via an exponential form, ie., 
that M{^,t) (X e"*^^"*"^^^-*. In fact this is not correct- the 
crucial conclusion in our paper about TdeiC) was that it 
scales as so that we can write M{^,t) as a function 



of one scaled variable alone, satisfying the scaling law 



Mi^,t)^M{t/e) = Miz), 



(1) 



where z = t/rde ^ i/C^- The actual form of M{z) is not 
so crucial. In our paper we used an exponential decay 
form for M{z) (this is just the standard r-approximation 
result, and so the simplest form to use). This was used 
for illustrative purposes only, and none of our conclusions 
is based on the specific form of this function- contrary 
to the comment's claim! In fact it is very easy to see 
that the functional form of M{z) is irrelevant for the 
final conclusion that the magnetization decay is ~ t^/'^, 
which simply and immediately follows from the use of the 
scaling form in M{t) = J d^M{^,t), ie., from the scaling 
form alone. This scaling form was found to be valid in our 
MC simulations for different lattice types, provided one 
assumes that t > Tq (ie., assuming initial transients are 
over) and that SM/M{t = 0) <C 1 (ie., assuming that the 
long-time multispin correlations have not yet set in). We 
also assumed that ^ £, Wd (here Wd is the half- 
width of the dipole energy distribution, is the nuclear 
spin bath parameter, and Tq is the characteristic single- 
molecule relaxation time^ hereafter we measure time in 
units of To and energies in units of ^o)- 

(ii) Alonso and Fernandez are basically arguing that 
the scaled variable z should be rather written as t/£}/P 
with the exponent p depending on the lattice type; cor- 
respondingly, M{t) ~ t^ and for FCC lattices they found 
p ~ 0.73. For SC lattices they found p = 0.5. In Fig.[Tlwe 
present MC solutions of kinetic equations for the magne- 
tization decay in the FCC lattice with 60^ spins (FCC, 
BCC and triclinic lattices were also analyzed in Ref.*). 
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After the usual initial transient behavior ~ t (already 
discussed in Ref.^ as well as^), we observe a long inter- 
val with the ^ ^/i relaxation at least until the fractional 
change in \M{t) - M(0)|/M(0) ~ 0.5. At longer times 
relaxation deviates from the ^ t^/^ behavior, but the du- 
ration of the ^ t^/^ interval in a demagnetized sample 
is longer than in a polarized sample (compare figures in 
Refs.2, andM). 
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FIG. 1: M{t) vs ^/i in the FCC lattice for initial magnetiza- 
tion Af (0) — 0.15 and three values of Ed' Ed ~ 2.5 (circles), 
Ed ~ 5 (diamonds) and Ed ~ 10 (triangles). The inset shows 
the effective relaxation time Tde((.) vs Solid lines in the 
inset are the shifted Lorentzian curves described in the text 
with a(ED) ~ 0.0415, 0.0205, 0.0114 for Ed = 2.5, 5, 10. The 
time-step in the MC simulations was St/ro — 10~^. 



Simultaneously, we find that the functional form of the 
effective relaxation rate in the FCC lattice is well de- 
scribed by the shifted Lorentzian T^^(f) = (l/7r)/(^^ + 
1) -I- a{ED) with small aiEo) ~ 1/Wd (again, = 1 hi 
our units); at small values of initial polarization in the 
FCC lattice Wd ~ 10i?D, where Ed is the strength of 
the dipole-dipole interactions between nearest neighbor 
spins. We did not find any qualitative difference between 
the SC, trichnic, BCC and FCC cases. 

It is important to notice that when the number of 
molecules in resonance is very small (as it is in the limit 
of very small values of £,o/Wd, particularly for FCC lat- 
tices), the usual initial transient of M(t) can be rather 
long, and its duration sensitive to the crystal structure 
because the dipole field spectrum is discrete in this limit. 
Then MC results on a small finite system are not rele- 
vant to experiments on macroscopic samples (where the 
number of molecules in resonance is also macroscopic). 

In the comment, and in a previous paper—, the authors 
have argued that they can fit the whole time range, in- 
cluding both very short and even infinitely long times, 
using a single theory with a single power law exponent p 
(see, eg.. Figs. 3,4 of the comment and Ref.®, paragraphs 
before Eqtn. 22 and after Eqtn. 24 of this paper). We 
find this implausible. Neither our MC simulationa^*^*^ 
nor our analytic work^ have ever claimed or attempted 
to explain more than the initial transient and the short- 
time relaxation that ensues after this transient, in the 
time interval before multi-spin correlations build up. It 
is in this restricted time range that we have argued for 
simple scaling and the associated square root relaxation. 
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